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ABSTRACT 

Several  general  formulas  relating  aerodynamic  forces  and  moments 
acting  on  finite  solid  bodies  immersed  in  a fluid  to  the  time-variation 
of  vorticity-moment  integrals  are  presented.  The  formulas  are  valid  for 
two-  and  three-dimensional  flows  and  are  shown  to  form  a theory  for  aero- 
dynamic forces  and  moments  which  encompasses  much  of  the  existing  aerodynamic 
theories . 


I.  INTRODUCTION 


The  problem  of  predicting  aerodynamic  forces  and  moments  acting  on 
finite  solid  bodies  immersed  in  and  moving  relative  to  a fluid  has  occupied 
the  center  stage  of  aerodynamic  research  from  the  end  of  the  nineteenth 
century  onwards.  In  fact,  it  is  this  focal  problem  that  distinguishes 
the  science  of  aerodynamics  from  other  branches  of  theoretical  fluid  mechanics. 
Studies  of  the  motion  of  the  fluid  relative  to  the  solid  bodies  of  course 
represent  a fundamental  aspect  of  aerodynamics.  In  the  majority  of  aerodynamic 
applications,  however,  such  studies  do  not  represent  ends  in  themselves. 

Rather,  they  are  undertaken  in  recognition  of  the  fact  that  the  motion 
of  the  fluid  is  ultimately  responsible  for  the  forces  and  moments  exerted 
on  the  solid  bodies  by  the  fluid. 

It  is  known  that  the  problem  of  finding  analytical  solutions  to  equations 
of  fluid  motion  associated  with  solid  configurations  of  practical  importance 
presents  considerable,  often  insurmountable,  mathematical  difficulties. 
Historically,  therefore,  the  most  remarkable  advances  in  aerodynamics 
were  brought  about  by  aerodynamic is ts  who  perceived  approaches  for  the 
prediction  of  aerodynamic  forces  and  moments  that  avoid,  as  much  as  possible, 
entanglement  with  the  details  of  the  fluid  motion.  For  example,  the 
Kutta-Joukowski  theorem,  i.e.  the  circulation  theory  for  two-dimensional 
steady  motion,  permits  the  lift  force  acting  on  a solid  body  to  be  determined 
from  a knowledge  of  the  circulation  about  the  body.  For  the  case  of  a 
two-dimensional  body  with  a sharp  trailing  edge,  the  Kutta  condition, 
requiring  the  .ear  stagnation  point  of  a potential  flow  to  be  located 
at  the  trailing  edge,  is  known  to  yield  acceptably  accurate  values  of 
the  circulation,  provided  that  the  flow  is  steady  and  does  not  separate 
over  an  appreciable  region  around  the  solid  body.  With  the  Kutta  condition. 
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the  problem  of  predicting  the  lift  force  for  such  a flow  is  reducible  to 
that  of  solving  an  integral  equation  (Ref.  1).  The  unknown  function  of  this 
integral  equation  is  a distribution  of  singularity  (sources,  sinks,  and  vorti- 
ces) over  the  body  surface.  There  ia  no  need  to  know  the  fluid  motion  away 
from  the  body  except  that  it  is  reasonably  represented  by  a potential  flow 
about  the  body. 

Extensions  of  the  circulation  theory  to  three-dimensional  and  to 
unsteady  flows  are  generally  based  on  the  concept  of  bound  vortices 
and  free  vortices.  The  bound  vortices  are  considered  to  move  with  the 
solid  body.  The  free  vortices  are  considered  to  be  shed  from  the  solid 
body  either  because  of  the  Helmholtz  theorem,  which  states  that  a vortex 
filament  cannot  begin  or  terminate  within  the  fluid  domain,  or  because  the 
total  circulation  of  the  entire  fluid  system  is  required  to  be  zero.  These 
free  vortices  are  located  at  finite  distances  from  the  solid  bodies  and 
they  make  a finite  contribution  to  the  fluid  motion  near  the  solid  surfaces. 
This  contribution  is  quantitatively  determinate  once  the  spatial  distribution 
of  the  free  vortices  is  known.  This  distribution  is  dependent  on  the  com- 
plex processes  of  shedding  of  vortices  and  subsequent  transport  of  free 
vortices  in  the  fluid.  To  avoid  the  detailed  computation  necessary  for  an 
accurate  determination  of  this  distribution,  most  previous  authors  have 
prescribed  the  motion  of  free  vortices  in  a simple  manner.  For  example, 
with  Prandti's  lifting-line  theory,  it  is  usually  assumed  that  the  free 
vortices  remain  stationary  relative  to  the  freestream. 

The  circulation  theory  and  its  extensions  permit  reasonably  accurate 
predictions  of  the  lift  force  for  some  solid  configurations  and  flow  environ- 
ments. The  scope  of  applicability  of  this  theory  and  its  extensions  has  not 
been  precisely  established.  It  is  well-known,  however,  that  considerable 
uncertainties  exist  regarding  the  value  of  the  circulation  to  be  prescribed 
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in  cases  where  the  solid  body  does  not  possess  a sharp  trailing  edge,  where 
appreciable  regions  of  flow  separation  exist,  or  where  the  motion  of  the  solid 
is  time-dependent.  These  uncertainties  arise  mainly  because  assumptions,  or 
hypotheses,  utilized  in  the  development  of  the  theory  are  often  ad-hoc  and, 
consequently,  the  theory  is  not  readily  interpreted  as  an  approximation  of 
a specific  physical  phenomenon.  For  example,  the  bound  vortex  is  usually 
described  in  well-known  treatises  as  "replacing"  an  airfoil  (or  a wing). 

The  vorticity,  defined  as  the  curl  of  the  velocity,  is  twice  the  angular 
velocity  in  a solid  region.  If  the  airfoil  is  not  rotating,  then  clearly 
it  does  not  possess  vorticity.  The  bound  vortex  is  therefore  not  an  approxi- 
mation of  the  airfoil.  To  the  experienced  aerodynamicist,  the  remarkable 
agreement  between  the  predicted  and  measured  lift  forces  under  certain  circum- 
stances points  to  the  inevitable  conclusion  that  the  circulation  theory  is  a 
reasonable  approximation  of  the  physics  of  the  problem  under  these  circum- 
stances. It  was  suggested  (Ref.  2)  that  the  bound  vortex  does  represent  the 


airfoil  plus  its  boundary  layers.  The  fact  that  such  an  interpretation  wus 
emphasized  so  very  recently  (in  1976)  has  motivated  the  principal  investigator 
to  examine  critically  the  historical  development  of  the  circulation  theory, 
together  with  its  more  recent  extensions,  with  the  hope  of  facilitating 
the  interpretation  of  the  circulation  theory  and  thus  contributing  to 
its  further  development  for  three-dimensional  flows,  unsteady  flows,  and 
separated  flows.  In  the  course  of  this  study,  several  general  formulas 
relating  aerodynamic  forces  and  moments  to  the  time-variation  of  vorticity- 
moment  integrals  were  uncovered.  It  became  evident  that  these  formulas  have 
far  reaching  consequences  in  the  realms  both  of  theoretical  aerodynamics  and  of 
computational  aerodynamics.  In  particular,  it  was  shown  that  the  formulas  form 
a theory  for  aerodynamic  forces  and  moments.  This  theory  encompasses  much  of 
the  existing  aerodynamic  theories.  For  example,  the  circulation  theory  for 
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steady  lift  and  its  various  extensions  are  readily  interpreted  as  various  levels 
of  approximation  of  the  general  theory.  The  purpose  of  this  report  is  to  present 
this  general  theory. 

A distinguishing  feature  of  the  present  theory  is  that  the  concept  of  bound 
vortex,  or  that  of  singularity  elements  such  as  sources,  sinks,  and  vortex  fila- 
ments, is  not  embodied  in  the  general  formulas  forming  the  theory.  Rather,  the 
actual  vorticity  distributions  of  the  fluid  and  of  the  solid,  the  latter  being  re- 
lated to  the  rotational  motion  of  the  solid,  enter  these  formulas.  The  starting 
point  of  the  present  theory  is  a rotational  flow  analysis.  Consequently,  the 
theory  is  applicable  to  viscous  flows.  This  freedom  from  "bondage"  is  important 
in  the  interpretation  of  the  various  aerodynamic  theories.  For  example,  it  permits 
a precise  definition  of  the  "circulation"  about  a two-dimensional  solid  in  the  case 
where  an  appreciable  region  of  separation  exists.  While  the  Kutta-Joukowski  theorem 
predicts  zero  drag,  the  general  formula  does  relate  the  time-variation  of  a vorti- 
city-moment  integral  to  a non-zero  drag,  including  the  profile  drag.  The  formulas 
predict  an  unsteady  drag  without  the  customary  energy  or  apparent  mass  consideration 
These  general  formulas,  in  fact,  clearly  point  out  the  basis  principles  for  mini- 
mizing the  drag  and  for  maximizing  the  lift.  Many  of  the  measures  proposed  for  drag 
reduction  and  for  lift-augmentation  are  readily  interpretable  on  the  basis  of 
these  principles.  The  present  research  deals  only  with  incompressible  flows, 
although  the  basic  principles  described  here  are  certainly  applicable  to 
compressible  flows  as  well. 

In  recent  years,  extensive  efforts  have  been  in  progress  at  many  research 
institutions  to  develop  numerical  methods  for  the  solution  of  aerodynamic  problems. 
An  ultimate  goal  of  these  efforts  is  to  make  available  methods  of  predicting 
aerodynamic  forces  and  moments  that  are  more  accurate  and  that  possess  a wider 
scope  of  validity  than  the  circulation  theory.  These  efforts  are  divisible  into 
two  major  categories.  In  one  category,  numerical  methods  are  being  developed 
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for  the  solution  of  inviscid  flow  equations.  The  free  vortices,  are  assumed  to 
convect,  but  not  to  diffuse,  with  the  fluid.  The  solid  body  is  represented  by 
a singularity  distribution.  Conceptually,  these  numerical  methods  utilize  the 
basic  assumptions  of  the  circulation  theory.  They  relax  the  restrictions  of 
analytical  methods  based  on  the  circulation  theory  through  detailed  compu- 
tation. For  example,  the  linearization  procedure  introduced  in  classical 
studies  of  the  unsteady  two-dimensional  airfoil  problem  (e.g.  Ref.  3)  is  no 
longer  necessary  if  numerical  methods  (e.g.  Ref.  4)  are  employed.  These 
numerical  methods  are  of  course  expected  to  be  subject  to  the  well-known 
limitations  of  the  inviscid  flow  assumption.  In  this  regard,  the  availa- 
bility of  the  general  formulas  are  expected  to  offer  clearer  interpretation 
of  these  numerical  methods  and  better  definition  of  their  scope  of  application. 

In  the  second  category,  numerical  methods  are  being  developed  for  the 
solution  of  differential  equations  governing  viscous  flows.  Impressive 
progress  has  been  made  in  recent  years  in  the  numerical  solution  of  two- 
dimensional  laminar  separated  flow  problems  as  well  as  in  the  establishment 
of  turbulence  models  for  separated  flows.  For  three-dimensional  separated 
flows,  the  development  of  numerical  methods  is  hindered  by  excessive  compu- 
tation requirements  (Ref.  5).  Methods  that  possess  superior  computational 
efficiency  are  therefore  of  critical  importance.  During  the  past  few 
years,  the  principal  ittvestigator  and  his  co-workers  have  developed  a 
new  numerical  approach  which  permits  the  confinement  of  the  solution  field 
to  the  vortical  region  of  the  flow  (Ref.  6).  In  comparison  to  available 
finite-difference  and  finite-element  methods,  all  of  which  require  the 
solution  field  to  include  the  potential  region  in  addition  to  the  vortical 
region,  the  new  approach  requires  a drastically  smaller  number  of  grid 
points.  The  computation  requirements  using  this  new  approach,  called  the 
integro-dif ferential  approach,  are  consequently  drastically  smaller  than 


those  using  other  methods  (Ref.  7,8).  The  new  approach  uses  the  vorticity 
vector  as  a field  variable  in  place  of  the  pressure.  The  general  formulas 
presented  in  this  report  provide  a convenient  means  of  computing  the  aero- 
dynamic forces  and  moments  directly  from  vorticity  distributions.  That  is, 
using  these  general  formulas,  it  is  no  longer  necessary  to  compute  first  the 
pressure  and  shear  stress  distributions  on  the  solid  surfaces  from  the  value 
and  normal  gradient  of  vorticity  on  the  surface  (Ref.  8)  and  then  the  inte- 
grated forces  and  moments.  In  addition,  the  general  formulas  suggested 
several  promising  techniques  for  minimizing  the  required  computation. 

Several  important  theorems  of  fluid  dynamics  are  utilized  in  the  deri- 
vation of  general  formulas  for  aerodynamic  forces  and  moments.  Similar 
theorems  and  formulas  are  given  in  many  well-known  textbooks  on  aerodynamics, 
e.g.  Refs.  9 and  10.  These  earlier  theorems  and  formulas,  however,  are 
traditionally  considered  in  the  context  of  an  infinite  limitless  fluid, 
i.e.,  an  infinite  fluid  with  no  internal  boundaries.  The  present  theorems 
and  formulas  are  valid  in  the  presence  of  internal  boundaries  that  repre- 
sent solid  surfaces.  The  practical  importance  of  this  more  general  validity 
needs  no  emphasis  since  the  interaction  between  the  solids  and  the  fluid  is 
indeed  what  the  subject  of  aerodynamics  is  about. 

The  subject  of  this  report  received  some  attention  in  the  1950's,  nearly 
a quarter  of  a century  ago.  Phillips  (Ref.  11)  presented  a formula  relating 
the  fluid  momentum  to  an  integral  of  vorticity  moment  for  a two-dimensional 
flow  associated  with  a cylinder  in  translation.  The  general  formulas 
presented  in  this  report  are  valid  for  two-  and  three-dimensional  flows 
associated  with  one  or  more  finite  solid  bodies  of  any  arbitrary  shape 
executing  any  prescribed  steady  or  time-dependent  translation  and/or  ro- 
tation. Moreau  (Ref.  12)  presented  formulas  for  a limitless  fluid  and 
for  a portion  of  fluid  subject  to  certain  "order"  conditions  at  infinity. 
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True 8 dell  (Ref.  13)  comnented  that  "Moreau  emphasized  his  application  to  a 
limitless  fluid,  all  but  a finite  interior  part  of  which  is  in  irrotational 
or  circulation  preserving  motion.  In  this  connection  we  should  beware  of 
the  extremely  strong  order  conditions  at  infinity  required  in  order  to 
get  simple  results,  order  conditions,  indeed,  which  possibly  may  never 
be  satisfied".  In  this  report,  formulas  are  rigorously  derived  for  the 
viscous  flow  of  fluids  past  finite  bodies,  using  order  conditions  at  in- 
finity that  are  shown  to  be  satisfied  under  quite  general  circumstances. 

The  applications  of  the  theory  formed  from  these  formulas  will  be  presented 
in  future  reports. 
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II.  VORTICITY  DYNAMICS 

The  time-dependent  motion  of  an  infinite  incompressible  fluid  with 
uniform  viscosity  relative  to  one  or  more  iranersed  solid  bodies  is  considered 
in  the  present  study.  The  solid  bodies  are  initially  at  rest  in  the  fluid 
also  at  rest  and  are  located  within  finite  distances  from  one  another. 
Subsequent  prescribed  motions  of  the  solid  bodies  induce  a corresponding 
motion  of  the  fluid.  At  large  time  levels  after  the  motion  has  initiated, 
if  the  solid  bodies  move  uniformly  at  a constant  translational  velocity 
relative  to  the  freestream,  then  the  possibility  of  an  asymptotic  steady 
flow  exists.  Alternatively,  the  possibility  of  a time-dependent  flow 
involving  periodic  vortex  shedding,  as  evidenced  by  the  well-known  Karman 
vortex  street  behind  a circular  cylinder,  also  exists.  In  the  present 
work,  a steady  flow,  when  it  exists,  is  considered  to  be  approached  asymp- 
totically at  large  time  levels  after  the  initiation  of  the  solid  motion. 

If  the  solid  bodies  do  not  move  uniformly,  or  if  the  solid  motion  is  time- 
dependent,  then  the  motion  of  the  fluid  is  necessarily  time-dependent. 

The  familiar  differential  equations  describing  the  time-dependent 
fluid  motion  are  the  continuity  and  Navier-Stokes  equations: 

V.v  = 0 (II-l) 

+ (v.  V ) v = - p »p+v  V v (II-2) 

where  v,  p,  p,  and  v are  respectively  the  velocity  vector,  the  pressure, 
the  density,  and  the  kinematic  viscosity  of  the  fluid,  p being  a constant 
in  the  present  study.  For  simplicity,  the  kinematic  viscosity  of  the 
fluid  is  considered  to  be  uniform  in  this  report.  It  is  not  difficult 


lam  if  d i f t n 


to  generalise  the  analyses  given  here  to  flows  where  the  viscosity  is 
not  uniform.  Such  a generalization,  however,  is  not  essential  to  the 
purpose  of  the  present  work. 

In  this  report,  the  region  occupied  by  the  fluid  is  designated  R^. 

A coordinate  system  with  its  origin  located  within  finite  distances  from 
all  solid  surfaces,  collectively  designated  by  Bg,  is  used.  Unless  otherwise 
specified,  this  coordinate  system  is  considered  to  be  at  rest  relative 
to  the  freestream.  The  fluid  region  is  bounded  internally  by  Bg 

and  externally  by  a close  boundary  ^ at  infinity.  The  region  occupied 
by  the  jth  solid  body  is  designated  Rj,  which  is  bounded  externally 
by  Bj.  The  limitless  region  jointly  occupied  by  all  the  solid  bodies 

and  the  fluid  is  designated  R . 

00 

It  is  convenient  to  introduce  the  vorticity  vector  u defined  by 


+ + -f 

V x v * b> 


(II-3) 


and  to  consider  the  vorticity  transport  equation 


■|^  « ^ X (yxZ  ) + V V1^ 


(II-4) 


obtained  by  taking  the  curl  of  both  sides  of  Eq.  (II-2)  and  using  Eqr. 

(II-l)  and  (II-3). 

The  set  of  equations  (II-l),  (II-3),  and  (II-4)  replaces  the  set 
of  equations  (II-l)  and  (IT-2).  There  are  several  reasons  for  using  ui 
in  the  formulation  of  the  problem.  In  the  first  place,  the  circulation 
theory  for  the  lift  force  suggests  that  the  vorticity  of  the  flow,  which 
ultimately  should  be  traced  to  the  circulation,  is  responsible  for  the 
forces  and  moments  exerted  by  the  fluid  on  the  solid.  Secondly,  as  shown 
in  Reference  6,  the  use  of  the  vorticity  vector,  which  is  intimately  connected 
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with  viscosity  effects,  permits  the  solution  field  for  the  incompressible 
flow  problem  to  be  confined  to  the  viscous  region  only.  Thirdly,  the 
set  of  equations  in  terar  of  u decomposes  conveniently  into  a kineaatic 
aspect  and  a kinetic  aspect,  each  aspect  constituting  an  entity  by  itself. 

The  first  feature  stated  above  provided  the  motivation  for  the  present 
effort  in  developing  general  formulas  relating  aerodynamic  forces  and 
moments  to  the  time-variation  of  vorticity-moment  integrals.  The  advantages 
offered  by  the  second  feature  have  been  extensively  studied  by  this  author 
and  his  coworkers  in  a series  of  previous  articles  in  the  context  of  computation 
methods.  Most  of  the  numerical  results  obtained  thus  far  have  been  for 
two-dimensional  incompressible  laminar  flows,  both  time-dependent  (Refs.  6,  7,  8) 
and  steady-state  (Ref.  14).  Some  results  have  been  obtained  recently, 
however,  for  relatively  simple  turbulent  flows  (Ref.  15).  In  addition, 
extensions  to  compressible  flow  problems  have  been  suggested  (Refs.  16,  17). 

The  importance  of  the  third  feature  in  the  present  work  is  due  to 
the  fact  that  the  physical  processes  of  flow  development  are  clearly  delineated 
once  the  overall  problem  is  decomposed  into  its  kinematic  and  kinetic 
aspects.  In  particular,  with  vorticity  as  a field  variable,  considerable 
insight  is  gained  by  examining  the  differential  equations  describing  these 
two  aspects,  without  employing  detailed  mathematical  or  mnerical  analyses. 

The  kinematic  aspect  of  the  problem  concerns  the  relationship  between 
the  vorticity  distribution  at  any  given  instant  of  time  and  the  velocity 
distribution  at  the  same  instant.  The  differential  equations  describing 
this  aspect  are  the  continuity  and  vorticity  definition  equations  (II -1) 
and  (£1-4).  Since  the  density  of  the  solid  bodies  does  not  undergo  appreciable 
change,  the  continuity  equation  (II-l)  ia  obviously  valid  within  the  solid 
regions  R.  as  well  as  in  the  fluid  region  R^.  If,  within  a solid  region 
R.,  one  defines  a vorticity  field  according  to  Eq.  (II-3),  then  the  kinematics 
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of  Che  solid  bodies  and  Che  fluid  are  described  by  Che  sane  differencial 
equaCions.  The  stress-strain  relations,  which  dif ferendaCe  che  fluid 
from  che  solid  bodies  kinetically,  do  noc  enCer  Che  kineutic  relacion 
beCweeen  Che  velociCy  field  and  Che  vorCiciCy  field.  As  a consequence, 

Che  solid  bodies  and  Che  fluid  can  be  CreaCed  together  as  one  kinenatical 
system. 

The  recognization  of  Che  fact  just  mentioned  makes  it  relatively 
simple  Co  derive  che  kinematic  theorems  and  formulas  presented  in  che 
report.  For  these  theorems  and  formulas,  Che  region  of  interest  is  limitless 
and  th«  differential  equations  leading  Co  these  theorems  and  formulas 
are  linear.  In  fact,  all  previously  available  kinematic  theorems  and 
formulas  derived  for  an  infinite  unlimited  fluid  are  immediately  applicable 
to  the  present  situation  of  an  infinite  fluid  with  one  or  more  immersed 
solid  bodies.  This  fact  is  not  well  recognized.  The  treatment  of  the 
solid  bodies  and  the  fluid  together  as  one  kinematical  system  has  not  received 
emphasis  in  the  literature.  In  this  report,  kinematic  theorems  and  formulas 
are  derived  by  treating  the  solid  bodies  and  the  fluid  together.  The 
author  has,  in  .addition,  re-derived  each  of  these  formulas  and  theorems 
by  considering  only  the  limited  fluid  region  which  is  bounded  internally 
by  solid  surfaces.  The  presence  of  boundaries  makes  the  derivations  lengthy 
and  algebraically  tedious.  In  this  report,  only  a few  of  these  re-derivations 
are  presented  to  demonstrate  the  validity  of  treating  the  solid  bodies 
and  the  fluid  together  as  one  kinematical  system. 

The  kinetic  aspect  of  the  problem  is  concerned  with  the  development 
of  the  vorticity  field  with  time.  This  aspect  is  described  by  the  vorticity 
transport  equation  (II-4).  This  equation  is  non-linear  in  the  sense  that 
the  first  term  on  its  right-hand  side  involves  the  product  of  v and 

■f  4 -*• 

ta>,  and  v is  kinematically  a function  of  u> . This  equation  is  valid 
only  in  the  fluid  domain,  which  is  limited. 
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Because  of  the  non-linearity  of  the  differential  equation  and  the 
necessity  of  treating  a 1 ini ted  region,  the  analysis  of  the  kinetic  aspect 
of  the  problem  presents  greater  smthematical  difficulties  than  the  analysis 
of  the  kinematic  aspect.  Undtr  certain  conditions,  it  is  possible  to 
specify  the,  vorticity  field  approximately  without  actually  solving  the 
vorticity  transport  equation.  It  is  then  only  necessary  to  deal  with 
the  kinenatic  aspect  of  the  problem  which,  as  stated  earlier,  is  described 
by  linear  differential  equations  applicable  to  the  entire  limitless  region 
^ . Under  these  circumstances,  the  aerodynamic  forces  and  moments  are 
sometimes  obtainable  in  a relatively  simple  manner. 

The  available  literature  on  classical  theories  for  aerodynamic  forces 
and  moments  shows  that  the  avoidance  of  the  kinetic  part  of  the  problem 
has  beet,  an  essential  ingredient  of  these  theories.  There  have  been  many 
recent  efforts,  via  numerical  methods,  to  treat  the  kinetic  aspect  either 
partially,  e.g.,  on  the  basis  of  the  inviscid  flow  equations,  or  fully. 

These  efforts  can  benefit  from  a clear  understanding  of  the  physical  processes 
involved  in  the  development  of  the  vorticity  field  in  the  fluid.  A significant 
amount  of  information  already  exists  in  the  literature  on  this  topic  (e.g. 

See  Ref.  10).  Those  features  of  vorticity-field  development  that  are 
pertinent  to  the  present  work  are  described  briefly  below. 

For  an  inviscid  fluid,  the  last  term  in  Eq.  (II-4)  vanishes  and  the 
vorticity  is  convected  with  the  fluid  in  the  sense  that  the  vorticity 
flux  u*ds  associated  with  each  material  element  ds  moving  with  the 
fluid  remains  a constant  for  all  times.  This  well-known  theorem  of  Relmholts, 
a proof  of  which  is  available  in  many  textbooks,  e.g.  Ref.  18,  is  modified 
in  the  case  of  a real  fluid  by  the  process  of  vorticity  diffusion.  According 
to  Eq.  (11-4),  changes  in  the  vorticity  flux  take  place  only  by  diffusion. 
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Vorticity  flux  cannot  be  created  or  destroyed  in  the  interior  of  a fluid. 

For  the  problem  under  consideration,  the  vorticity  is  obviously  everywhere 
zero  prior  to  the  impulsive  start  of  the  siotion  of  the  solid  bodies. 

The  interior  of  the  fluid  domain  therefore  can  become  vortical  only  if 
vorticity  diffuses  across  the  boundaries  of  the  fluid  region.  Consequently, 
inaediately  after  the  onset  of  the  motion,  the  vorticity  is  everywhere 
zero  in  the  fluid  except  at  the  boundaries  in  contact  with  the  solid  bodies. 
That  is,  the  fluid  motion  immediately  after  the  onset  of  the  motion  has 
a non-zero  tangential  velocity  relative  to  the  solid  bodies  at  the  solid 
boundaries.  The  discontinuity  in  tangential  velocity  constitutes  a sheet 
of  concentrated  vorticity  (vortex  sheet)  at  the  boundaries.  At  subsequent 
time  levels,  this  concentrated  vorticity  spreads  into  the  interior  of 
the  fluid  domain  by  diffusion  and,  once  there,  is  transported  away  from 
the  boundaries  by  both  convection  and  diffusion.  At  the  same  time,  the 
no-slip  condition  provides  a mechanism  for  the  continual  generation  of 
vorticity  at  the  boundaries.  The  general  flow  pattern  therefore  contains 
vortical  regions  surrounding  the  solid  bodies  and  vortical  wakes  trailing 
the  solid  bodies.  Outside  of  these  vortical  regions  and  wakes,  the  flow 
is  essentially  free  of  vorticity  and  therefore  irrotational . In  particular, 
if  the  flow  Reynolds  number  is  not  small,  then  the  vorticity  spreads  by 
diffusion  only  a short  distance  from  the  boundaries  before  it  is  carried 
away  with  the  fluid  by  convection.  Therefore  a large  region  of  the  fluid, 
ahead  and  to  the  side  of  the  solid  bodies,  is  essentially  free  of  vorticity 
and  irrotational. 

In  the  next  two  Chapters,  a number  of  theorems  and  formulas  are  derived 
using  the  above  described  kinematic  and  kinetic  characteristics  of  the 
flow. 
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III.  SELECTED  THEOREMS  AND  FORMULAS  FOR  THE  KINETIC  ASPECT 


1.  Asymptotic  Behavior  and  Effective  Extent  of  Vorticity  Field. 

For  the  present  problem,  the  vorticity  decays  exponentially 

j 

with  increasing  distance  from  the  origin  at  large  distances  for 
all  finite  time  levels  after  the  onset  of  the  motion  and  the  effective 
extend  of  the  vortical  region  is  finite.  The  finite  extent  of  the 
vortical  region  is  a consequence  of  the  fact  that  the  vorticity 
is  transported  in  the  fluid  by  finite-rate  processes  and  cannot 
be  created  in  the  interior  of  a fluid  domain.  This  fact  as  well 
as  the  asymptotic  behavior  of  the  vorticity  distribution  at  large 
distances  from  solid  bodies  are  established  below  by  considering 
the  fundamental  solution  of  the  diffusion  equation. 

The  fundamental  solution  F of  the  diffusion  equation,  i.e.  the 
Green's  function  for  an  infinite  unlimited  region,  can  be  expressed  as 
(Ref.  18  and  19) 


F(r, 


1 

d 

(4HV(t-to))2 


exp 


♦ ♦ 2 
/ I r"ro| 

1 4 v?t=t“y 

o 


(III-l) 


where  r and  tq  are  position  vectors  and  d is  the  dimensionality 

of  the  problem,  i.e.,  d * 1,  2,  and  3 respectively  for  one-,  two-,  and 

three-dimensional  problems.  For  the  present  problem,  the  fundamental 

solution  represents  the  vorticity  distribution  at  the  time  level  t resulting 

from  the  diffusion  of  a concentrated  vorticity  which  is  of  "unit"  strength 

and  is  located  at  the  point  ^ at  the  time  level  t . with  t <t. 

r o o o 


15 


If,  at  the  tine  level  tQ,  the  vortieity  is  non-zero,  only  in  an  elemental 

region  dR  located  at  r and  the  value  of  vortieity  in  dR  is  u . 

o o o o 

-► 

then  the  vortieity  distribution  at  the  subsequent  time  level  t is  F*>odRQ. 
If,  at  the  time  level  tQ,  the  vortieity  distribution  is  known  in  the 
unlimited  region  R^,  then  the  vortieity  distribution  at  a subsequent 
time  level  t is  expressible  as  an  integral: 


ui  (r , t)  ■ / F*Ki 


t )dR 
o o 


(III-2) 


where  the  subscript  for  dRQ  indicates  that  the  integration  is  performed 

in  the  r space, 
o 

Equation  (III-2)  is  valid  in  an  infinite  unlimited  region  in  the 
absence  of  convection.  The  form  of  the  fundamental  solution  F clearly 
shows  that,  if  the  vortieity  is  non-zero  at  the  time  level  tQ  only  within 
finite  distances  from  the  origin,  then  the  vortieity  at  any  subsequent 
finite  time  level  t approaches  zero  exponentially  with  increasing  distance 
r from  the  origin,  at  large  distances.  Therefore,  the  vortical  region 
is  effectively  confined  to  a finite  region  at  any  finite  time  level  t. 

For  a fluid  region  bounded  internally  by  solid  surfaces  and  in  which 
convective  process  is  present,  Eq.(lII-2)  needs  to  be  generalized.  It 
is  obvious  that  the  convective  process,  being  one  of  finite  rate,  does 
not  alter  the  above  conclusions  regarding  the  effective  extent  of  vortical 
regions  and  the  asymptotic  behavior  of  vortieity  at  any  finite  time  level. 
Similarly,  the  presence  of  solid  surfaces  in  the  fluid  provides  a mechanism 
for  introducing  vortieity  at  the  boundaries  of  the  fluid  region  and,  as 
long  as  these  boundaries  are  within  finite  distances  from  the  origin, 
the  introduction  of  vortieity  does  not  alter  the  above  conclusions.  For 
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two-dimensional  flows,  the  generalised  version  of  Eq.  (1H-2.)  is  expressible 
as  (Refs.  17  and  20): 


where  the  subscript  "o"  indicates  that  the  variables,  differentiations, 

♦ + 

and  integrations  are  in  the  rQ,  tQ  space,  e.g.  ^ ■ <i)0'r0»  The 

second  and  third  integrals  of  Eq.  (III-3)  represent  respectively  the  convective 

process  and  the  effect  of  solid  boundaries.  Each  of  the  integrands  in 

-*■ 

the  integrals  of  Eq.  (III-3)  is  directly  proportional  to  P and/or  V F. 

At  any  given  time  level  t and  at  large  distances  r from  the  origin, 

F and  $ F decays  exponentially  with  increasing  r.  It  follows  that, 
for  any  problem  in  which  Bg  consists  of  finite  surfaces  located  within 
finite  distances  from  the  origin,  if  the  vorticity  decays  exponential 
with  increasing  r at  large  r for  all  time  levels  previous  to  t,  then 
the  vorticity  decays  exponentially  with  increasing  r also  at  the  time 
level  t.  In  fact,  for  such  a problem,  if  the  vorticity  is  confined  within 
finite  distances  from  the  origin  at  any  given  instant  of  time,  then  the 
vorticity  decays  exponentially  with  increasing  r at  large  r at  all 
subsequent  finite  time  levels.  For  the  present  problem  involving  an  impulsively 
started  motion,  the  vorticity  is  non-zero  immediately  after  the  onset 
of  the  motion  only  on  the  solid  surfaces.  Consequently,  the  statements 
made  in  the  first  paragraph  of  this  Section  are  true  for  two-dimensional 
flows.  It  can  be  similarly  shown  that  these  statements  are  also  true 
for  three-dimensional  flows. 
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2.  Principle  of  Total  Vorticity  Conservation  — Two-Dimensional  Flows. 

It  shall  be  shown  that,  for  the  flow  of  an  incompressible  fluid  past 
solid  bodies,  the  total  vorticity  in  the  infinite  unlimited  space  occupied 
jointly  by  the  fluid  and  the  solid  bodies  is  invariant  with  respect  to 
time,  provided  that  an  order  condition  for  the  vorticity  at  infinity  is 
satisfied.  This  order  condition  is  that  the  vorticity  approaches  eero 
as  r~°,  where  n>d,  d being  the  dimensionality  of  the  problem. 

The  above  statement  will  be  referred  to  as  the  principle  of  total 
vorticity  conservation  and  is  expressible  mathematically  as 


0 


(III-4) 


provided  that  u>  approaches  r n as  r >•,  with  n>d. 

As  noted  earlier,  the  solid  bodies  and  the  fluid  can  be  treated  together 
as  one  kinematic  system  in  the  present  problem.  For  three-dimensional 
f lows , since  the  vorticity  field  is  solenoidal  and  is  effectively  confined 
to  a finite  region,  all  vorticity-lines  in  the  combined  system  form  closed 
curves.  Consequently,  one  has 


l 


(i)  dR 


(III-5) 


In  the  next  Chapter,  a proof  of  Eq.  (HI-5)  is  presented  treating  the 
solid  bodies  and  the  fluid  as  separate  kinematicai  systems.  Clearly,  Eq.  (Ill-A)j 
follows  directly  from  Eq.  (III-5).  Thus  the  total  vorticity  is  not  only 
conserved,  it  must  always  be  eero  in  three-dimensional  flows. 

Por  two-dimensional  flows,  the  vorticity-lines  are  directed  perpendicular 
to  the  plane  of  the  flow.  The  vorticity-lines  extend  to  infinity  in  the 
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direction  perpendicular  to  the  flow,  and  they  do  not  form  closed  curves 
in  the  plane  of  the  flow.  The  total  vorticity  of  the  fluid  and  the  solid 
bodies,  that  is,  the  integral  of  vorticity  in  the  infinite  unlimited  plane 
of  the  flow,  is  still  conserved.  This  principle  of  total  vorticity  conservation 
for  two-diqpensional  flows  is  not  usually  discussed  in  standard  treatises 
on  fluid  dynamics.  A proof  of  this  principle  for  two-dimensional  flows 
is  given  below  on  the  basis  of  the  kinetics  of  the  problem. 

For  two-dimensional  flows,  the  vorticity  transport  equation  (11-4) 
can  be  rewritten  as 


— * -vV  *7xu 


(II 1-6) 


where  ^ denotes  a substantial  derivative. 

The  time  rate  of  change  of  the  total  vorticity  of  the  fluid  is 


(III-7) 


where  R^(t)  is  the  entire  region  occupied  by  the  fluid  and  is  a function 
of  the  time. 

Placing  Eq.  ( I II— 6 ) into  Eq.  (III-7)  and  using  Stoke 's  theorem, 
one  obtains 
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(III-8) 


where  it  is  a unit  outward  directed  normal  vector. 

In  Eq.  (1II-8),  the  boundary  B consists  of  the  . u-.d  boundary 
Ba  and  the  boundary  at  infinity  B^.  The  contribution  ci  Ba  to  Eq.dll  5) 
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is  zero  provided  that  y approach  zero  as  r for  large  r,  with  o >2. 
Since  u approaches  zero  exponentially  with  increasing  r for  large 
r,  the  above  order  condition  at  large  r is  satisfied  for  the  present 
problem.  The  boundary  B in  Eq.  (III-8)  can  therefore  be  replaced  by 
the  solid  boundary  Bg. 

Equation  (II-2)  can  be  rewritten  as 


. -I^p-vVxu 
Dt  p 

Taking  the  vector  product  of  each  term  in  Eq.  (III-9)  with 
integrate  the  resulting  equation  around  Bg,  one  has 


(III-9) 
n and 
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The  first  integral  on  the  right-hand  side  of  Eq.  (111-10)  is  zero 
by  virtue  of  the  single-valuedness  of  pressure  on  Bg.  This  fact  can 
also  be  shown  by  using  ..he  Stoke 's  theorem  and  the  fact  that  the  curl 
of  the  gradient  of  any  scalar  function  is  zero.  Combining  Eqs.  (III-8) 
and  (III-10)  therefore  gives 


fc/  ■«  - - /S ***•  (m-ii) 
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Consider  now  the  region  Rc  bound  externally  by  B . With  the  no- 
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slip  condition,  the  substantial  acceleration  ^ on  Bg  is  identical 
for  the  solid  bodies  and  for  the  fluid.  Using  Stokes's  theorem  and  the 
fact  that  the  outward  normal  vector  for  Rg  is  directed  opposite  to  that 
for  Rf,  one  obtains: 

^ ; dR  - - / $ x(j£)  dR  (III-12) 

Rf(t)  Rg(t) 
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In  the  solid  region  Rj , the  velocity  vector  is  given  by 


v » v.  + W . x r 
J 3 


(III-13) 


where  v.  is  the  rectilinear  velocity  of  the  solid  body  j and  2L 
is  its  rotational  velocity.  For  two-dimensional  problems,  ^ j is  directed 
perpendicular  to  the  plane  of  the  flow.  The  vorticity  in  the  solid  body 
j is  readily  obtainable  by  taking  the  curl  of  the  velocity  vector  as 
given  by  Eq.  (111-13)  and  is 
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Using  vector  differential  identities,  it  is  simple  to  show  that 


in  each  of  the  solid  regions  One  thus  obtains  from 


Eq.  (III-12) 
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Equation  (III-15)  is  equivalent  to  Eq.  (III-4).  It  states  that  the 
rate  of  change  of  the  total  vorticity  is  equal  in  magnitude  and  opposite 
in  sign  to  the  rate  of  change  of  the  total  vorticity  in  the  solid  bodies, 
or,  equivalently,  the  total  vorticity  in  the  infinite  unlimited  region 
is  zero.  If  the  rotational  velocities  of  the  solid  bodies  are 

prescribed  functions  of  time,  then  the  rate  of  change  of  the  total  vorticity 
in  the  fluid  can  be  calculated  using  the  following  simple  formula: 


J u>  dR  ■ 
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where  N is  the  total  number  of  solid  bodies  present,  and  Rj  is  the 
size  of  the  solid  body  j. 

Equation  (III-15)  can  be  integrated  with  respect  to  time,  yielding 


r -*  r ■* 

J u>dR+J  0)  dR  • A 

Rf(t)  Rg(t)  (111-17) 


where  A is  a constant  vector. 

For  a motion  starting  from  rest,  the  total  vorticity  in  the  combined 
fluid  and  solid  regions  is  zero  before  the  onset  of  the  motion.  Consequently, 
A * 0 and  one  obtains  (II1-5),  which  states  that  the  vorticity  in  the 
combined  solid  and  fluid  regions,  i.e.,  the  infinite  unlimited  region, 
is  always  zero.  Thus,  if  the  solid  motion  is  prescribed  at  any  tisw  level, 
the  total  vorticity  in  the  fluid  is  easily  calculated  from 

+ N + 

J w dR  - -2  jSj  fl  jRj  (III-18) 
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There  are  several  conceptual  differences  between  the  principle  of 
total  vorticity  conservation  discussed  here  and  the  usual  understanding 
of  invariance  of  vorticity  integral  with  respect  to  time. 

In  the  present  work,  the  solid  regions  are  included  in  the  evaluation 
of  the  total  vorticity.  The  integrands  of  Eqs.  (III-4)  and  (III-5)  are 
piecewise  continuous  and  the  integrals  converge.  The  meaning  of  the  total 
vorticity  in  the  unlimited  infinite  space  jointly  occupied  by  the  fluid 
and  the  solid  bodies  is  unambiguous.  While  previous  discussions  of  the 
invariance  of  the  total  vorticity  are  usually  made  in  the  context  of  an 
inviscid  fluid  or  of  an  unlimited  infinite  fluid  region  (in  the  absence 
of  internal  boundaries),  the  present  study  utilizes  the  no-slip  condition 
at  the  solid  boundaries  and  permits  the  presence  of  such  boundaries  in 
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the  fluid.  In  Ref.  21,  Section  B.2,  a formula  similar  to  Eq.  (111-18), 
but  specialized  to  a single  solid  and  including  an  additional  term  repre- 
senting the  contribution  of  the  velocity  at  a surface  enclosing  the  solid,  is 
presented.  In  the  present  work,  it  is  shown  that  this  contribution  is 
absent  if  the  surface  is  sufficiently  distant  from  the  solid  bodies  and 
the  order  condition  for  the  vorticity  vector  is  satisfied.  The  formula 
(III-18)  is  derived  for  one  or  more  solid  bodies  in  the  present  Chapter. 

Ref.  21  also  gives  a formula  for  the  rate  of  change  of  the  total 
vorticity  in  the  fluid.  That  formula  contains  a term  representing  "conduction 
of  vorticity  through  the  solid  boundary."  It  is  pointed  out  that  further 
dynamic  (kinetic)  equations  are  needed  to  evaluate  this  term.  The  present 
result,  Eq.  (III-16) , shows  that  this  term  is  given  simply  by  the  rate 
of  change  of  the  total  vorticity  of  the  solid  bodies. 

In  Ref.  10,  it  is  shown  that  in  three-dimensions  the  total  vorticity 
is  zero  in  a region  which  contains  the  fluid  region  and  "a  region  extending 
beyond  the  actual  boundaries".  The  present  results  show  that  the  proper 
extension  of  the  fluid  region  is  simply  the  solid  regions  in  which  the 
correct  vorticity  values  to  assign  are  the  actual  vorticity  of  the  solid 
bodies.  For  two-dimensional  flows,  the  literature  emphasizes  the  possibility 
of  the  existence  of  a non-zero  circulation  around  closed  paths  at  large 
distances  from  the  solids.  According  to  the  present  results,  this  possibility 
does  not  exist  for  a real  fluid  at  any  finite  time  level  after  a motion 
has  started  from  rest.  An  evaluation  of  the  conduction  of  vorticity  through 
solid  boundaries  is  not  necessary  since  this  "conduction  process"  conserves 
the  total  vorticity  in  the  infinite  unlimited  region  occupied  jointly 
by  the  fluid  and  the  solid  bodies. 
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3.  Stress  Outside  the  Vortical  Regions. 

Outside  the  vortical  regions,  the  vorticity  is  zero  and  the  viscous 
stress  is  absent.  The  momentum  equation  (II-2)  simplifies  to 


P - - V(p  + p v2/2) 


(111-19) 


The  absence  of  vorticity  implies  the  existence  of  a scalar  potential 
9 such  that 


v » - V4>  (III-20) 

Placing  Bq.  (111-20)  into  Eq.  (111-19)  and  integrating  the  resulting 
equation  in  space  gives 


2 

P * p -§£  - 2%  + f(t)  (III-21) 

This  well-known  equation  for  unsteady  inviscid  pressure  will  be  utilized 
to  derive  general  formulas  for  aerodynamic  forces  and  moments  in  Chapter  V. 

It  should  be  noted  that  the  scalar  potential  is  single- valued.  For  three- 
dimensional  flows  past  finite  solid  bodies,  the  region  in  which  9 exists 
is  singly-connected  and  therefore  $ is  independent  of  path.  For  two- 
dimensional  flows,  the  region  in  which  $ exists  is  multiply-connected. 
However,  the  cyclic  constant  for  9 is  zero  since  the  vortical  regions 
are  of  finite  extent  and  the  total  vorticity  is  zero. 


] 
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IV.  SELECTED  THEOREMS  AND  FORMULAS  FOR  THE  KINEMATIC  ASPECT 


The  proof?  of  the  theorems  and  formulas  for  the  kinematic  aspect 
of  the  external  flow  problem  utilise  only  the  kinematic  relationship  between 
the  velocity  and  the  vorticity  fields  and  the  order  conditions  for  the 
vorticity  field. 

1.  Principle  of  Total  Vorticity  Conservation — Three-Dimensional  Flows. 


Let  n,  bj,  and  b2  be  a right-handed  set  of  orthogonal  vectors  on 
the  boundary  Bg.  bj,  b2  are  tangential  unit  vectors.  Let  the  velocity 
components  in  the  n,  b^,  and  b2  directions  be  vn>  v^j,  and  v^2  respectively. 
With  the  no-slip  condition,  these  velocity  components  are  identical  on 


Bg  for  the  fluid  and  for  the  solid  bodies.  The  normal  component  of  the 


vorticity  vector  on  the  boundary  Bg  is  given  by  b2  - bl  and 


3b  i 3b  j 


is 


continuous  across  Bg.  That  is  u)  . n on  Bg  for  the  fluid  is  identical 


to  that  for  the  solid  body. 

-*• 

Using  the  fact  that  ui  is  solenoidal,  one  obtains  through  vector 
differencial  identities: 


V • (.(*•*)  <*»  J * a.  u 


CIV-1) 


where  a is  an  arbitrary  constant  vector. 
The  divergence  theorem  then  gives 


J co  dR  « fi  w»n)dB 
R B 


(IV-2) 


Let  R be  the  fluid  region  R{.  The  boundary  B then  consists 
of  the  aolid  boundary  Bg  and  a boundary  at  infinity.  The  contribution 
of  the  boundary  at  infinity  ia  aero  provided  that  u>  approachea  aero 
aa  r n for  large  r,  with  n>3.  This  order  condition  ia  certainly 
net  in  the  present  problem  for  which  w decaya  exponentially  with  r 
for  large  r.  Therefore,  the  boundary  B in  Eq.  (IV-2)  can  be  replaced 
by  Bg. 


Since  the  normal  component  of  u ia  continuous  acroaa  Bg,  one  may 
replace  u>  . n in  the  right-hand  aide  of  Eq.  (IV-2)  by  2 H^.n,  where 
W.  ia  the  angular  velocity  of  the  j aolid  body.  One  then  haa 

J 


-*■  N +*■  + 
JudR  - 2 I ^ r(  fiy  n)dB 


Rf  J ‘ B, 

* 4 

Consider  now  the  region  R in  Eq.  (III-2)  to  be  the  region  Rj 
bounded  externally  by  Bj.  One  haa 


(IV-3) 


2 J r(  fl..n)dB  - 2 J $ . 
n.  J o J 


dR 


(IV-4) 


'j  “j 

Placing  Eq.  (IV-4)  into  Eq.  (IV-3)  and  noting  that  the  normal  vector  n in 
Eq.  (IV-3)  is  directed  outward  from  Rf  while  it  ia  directed  inward  from  Rf  in  Eq. 
(IV-4),  one  obtaina 


f w dR  - -2  I fi.R. 
J 4 J J 

R*  J 


(IV-5) 


which  ia  equivalent  to  Eq.  (III-5). 

2.  Biot-Savart's  Law. 

It  ia  ahown  in  Ref.  6 that,  by  uaing  the  fundamental  aolution  of 
the  Poiaaon'a  equation,  the  kinematica  of  the  problem,  i>e.  Eqa.  (II-l) 
and  (11-3), ia  expresaible  in  the  form  of  an  integral  representation: 


where  A and  d are  constant*  depending  on  the  dimensionality  of  the 
problem  A * An  and  d ■ 3 for  three-dimensional  problems,  A * 2n 

and  d • 2 for  two-dimensional  problems and  the  subscript  "o"  indicates 

that  the  variables  and  the  integrations  are  in  the  xq  space,  i.e.  w,  ■ w(rQ,  t) 
etc. 

Let  the  region  R be  the  fluid  region  R^  bounded  internally  by 

the  solid  boundary  B#  and  externally  by  a boundary  at  infinity.  With 

a coordinate  system  attached  to  the  freestream,  the  velocity  at  infinity 

is  *ero  and  the  contribution  of  the  boundary  integral  to  the  velocity 

field  is  therefore  limited  to  the  solid  boundary  B_. 

s 

The  contribution  of  the  solid  boundary  Bj  to  the  velocity  field 
is  expressible  as 

" $ ((VV^P  V dBo  C IV— 7 ) 

Bj 

where  P is  the  fundamental  solution  of  the  Poisson's  equation  and  is 
defined  by 


1 


three-dimensional  problems 
two-dimensional  problems 


Using  vector  identities,  Eq.  (IV-7)  can  be  rewritten  as 
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1.  ■ $ f Pv’  ,n  dB.  x f Pv  x n dB 

1 O O O _ • o o i 


o o o 


(IV-8) 


Consider  now  the  solid  region  Rj  bounded  externally  by  Bj.  Using 
the  divergence  theorem  and  Stoke 'a  theorem,  and  noting  that  the  outward 
normal  for  the  region  Rj  is  directed  opposite  to  n,  one  obtains  from 
Eq.  (IV-8) 


/ M(v(p  v)+  * x?o  x(P  V ) dRc 


(IV-9) 


Using  vector  differential  identities,  the  integrand  in  Eq.  (IV-9) 
can  be  re-expressed  as 


- 0 * vQ)  x?  P - vQ(  $QP)  (IV-10) 


It  is  easy  to  show  that 


y . $QP  • - V 2p  0 for  r£ 


Thus,  in  the  fluid  region  R_,  one  has, 


I.  - - 
3 


fij x(?-r) 


r 

| ■»  4 

RJ  lro"r 
3 


-a dR 

,d  o 


One  therefore  obtains  from  Eq.  (1V-6), 


x (ro-r) 


vV.t)  - -I  fiS 

L ivr 


3 o 


(IV-ll) 


I 


LJ. 


.<■  -*  • „v  f * * 


V ^ < * r 
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where  R is  Che  unlimited  infinite  region  jointly  occupied  by  the  fluid 
and  the  aolid  bodies.  Equation  (IV- 11)  is  an  expression  of  the  Biot-Savart's 
law  for  a distribution  of  vorticity  in  an  unlimited  infinite  region. 

In  fact,  by  considering  the  solid  bodies  and  the  fluid  together  as  one 
kinemacical' system,  Eq.  (IV-  11)  can  be  inmediately  written  down  without 
going  through  the  intermediate  steps  dealing  with  the  boundary  conditions 
on  Bs. 

3.  Asymptotic  Behavior  of  Velocity  Field. 

To  examine  the  asymptotic  behavior  of  the  velocity  field  at  large 
distances  from  the  solid  bodies,  one  re-expresses  Eq.  (IV-11)  as 


v(r,  t)  - V x f w pdR 
o o 


'l 


(IV-12) 


For  three-dimensional  problems,  expressing  the  function 
a Taylor  series  about  the  point  tq  ■ 0,  one  obtains 


I Vr  I 


as 


1 


i^nlrV  r*°-[<5°TCMr-o] 


Or 


+ i [(  V V ‘Y  ~)  I ”]  + 

?Vo:.  L ° °|VrflrJ 


■ ~ - ro.  + terms  of  order  r n,  n^  3 


I vr  I 


(IV-13) 


One  has  from  Eqs.  (IV-8),  (IV-12),  and  (IV-13),  the  following  expression 
for  the  velocity  vector 


<-  ■ '■»»  * « ■ v.  •>, . 


* 


' 
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i(r,  t)  - ^i).r  v»0  - y>  •/(»;.»  i>  v«, 

Rflo  » oo. 


(IV- 14) 


• + terms  of  order  r , 4 

The  first  term  on  the  right-hand  side  of  Eq.  (IV-4)  is  zero  because 
the  total  vorticity  in  the  infinite  unlimited  region  is  zero,  i.e.  Eq.(lV-5) 
or  Eq.  (III-5). 

Using  the  identity 

V ( (r  $ :/)  * ($  7<Sn)  r * (^7*  r ) w (IV-15) 

o v o r o o r o o r o o 


the  second  term  on  the  right-hand  side  of  Eq.  (IV-14)  is  rewritten,  with 
the  help  of  the  divergence  theorem,  as 


(IV-16) 


■Uli  . r )S  - (Vi  . Sfr  ) 

4 it  ( 2 J r o wo  r oo 


- j f (r  . V -)  r («»)  n )dB 
i J o r ooo  o 


The  last  integral  in  the  above  expression  vanishes  if  o>  approaches  zero  as 
r n for  large  r,  with  n>4.  This  condition  is  certainly  met,  as  discussed 
in  Chapter  III.  The  remainder  of  the  above  expression  can  be  rewritten 


K^7*«o> 


(IV-17) 


where  aQ  is  the  total  first  moment  of  the  vorticity  oj  defined  by 
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■*  <*♦ 

a - ( r x w dR 

o Jo  “o  o 


(IV-18) 


Equation  (IV-14)  is  thus  expressible  as 


v^»  ^ 3 ($  ^ . aQ)  + terms  of  order  r-n,  A (IV-19) 

The  first  term  on  the  righthand  side  of  Eq.  (IV-18)  is  of  order  r“3. 
Provided  that  is  not  zero,  v is  of  order  r~3  for  large  r. 

For  two-dimensional  problems,  the  function  In  can  be  expressed 

m , Ir0r! 

as  a Taylor  series  about  the  point  r * 0: 


In 


,r  -ri 

I ° ! 


i ? i 


In  — - rQ.  V In  — + terms  of  order  r~  , n>3 


Following  similar  procedures  as  that  given  above  for  three-dimensional 
problems,  one  obtains  for  two-dimensional  problems: 


v(r,  t) 


k (M> x I “odRo  - r **/<v  *ln  ?>  kdRc 

^ __  K ar>- 


+ terms  of  order  r~n,  n>3 


(IV- 20) 


The  first  term  on  the  right-hand  side  of  Eq.  (IV-20)  vanishes  because 
of  Eq.  (III-18).  The  second  term  may  be  rewritten,  by  noting  that 


4(V$ln7>  ~ ro(  2o.  $ In  i) 


- (7  In  -^)  x (r  x u ) 
r o o 


and  that  in  two-dimensional  problems 


The  first  term  on  the  righbhand  side  of  the  above  expression  vanishes 


for  r * 0.  One  thus  has,  for  two-dimensional  problems 


Therefore,  provided  that  a 


The  above  conclusions  for  two-dimensional  and  three-dimensional  problems 


can  be  combined  into  the  statement  that,  provided  that  - ® is  not  xero 


for  large  r,  d being  the  dimensionality  of  the 


problem 


4.  Velocity  Integrals  in  Large  Bounded  Regions 


Using  Eqs . (II-l)  and  (11-3)  and  vector  differential  identities 


one  has,  for  three-dimensional  problems 


The  divergence  theorem  and  its  corollaries  then  give 


J(r  x u>-2v)dR 


The  above  equation  can  be  re-expressed  as 
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Let  R be  the  region  Rj  occupied  by  the  solid  body  j.  Equation 
(IV-23)  gives 


1 r •* 

2 ? r X 

B. 

J 


(n  x v)dB 


- J 

R. 

3 


(IV-24) 


Consider  now  R in  Eq.  (IV-23)  to  be  R'^,  the  part  of  fluid  region 
bounded  internally  by  the  solid  surface  Bg  and  externally  by  a large 
but  finite  spherical  surface  r = L.  Let  this  spherical  surface  be  designated 
B . Let  L be  sufficiently  large  so  that  Bt  encloses  all  the  solid 
bodies  present.  The  boundary  B in  Eq.  ( IV-23)  now  consists  of  the  solid 
boundary  Bg  and  the  spherical  boundary  B^.  With  the  no-slip  condition, 
the  contribution  of  Bg  to  Eq.  (IV-23)  is  given  by  Eq.  (IV-24).  One 

■f 

then  has,  noting  that  for  R'^  the  normal  vector  n is  directed  into 
the  solid  whereas  in  Eq.  (IV-24)  n is  directed  into  the  fluid. 


"vdR 


-■  f,  / v.dR  * ~ J Cr  xw  )dR  ♦ 

J 1 R.  3 2 R’ 

J f 


N 


f r .dR 

R 3 

3 


-*■ 

f r x(n  x v)dB 


(IV-25) 


Equation  ( IV-25)  can  be  re-expressed  as 


xu>  )dR  - ~ 


x(n  x v)dB 


(IV-26) 
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where  is  the  spherical  region  r < L and  includes  the  fluid  region 

Rf'  as  well  as  all  the  solid  regions  Rj . 

Let  L be  sufficiently  large  so  that  the  total  vorticity  outside 
is  negligible  and  the  velocity  on  *L  ” accurately  given  by  the  firat 
term  on  the  right-hand  side  of  Eq.  (IV-19).  The  integral  over  the  surface 

in  Eq.  (IV-26)  is  finite  and  expressible  in  teraw  of  the  vorticity-monent 
integral  a as  shown  below. 

Let  i,  j,  and  k be  a set  of  right-handed  normal  unit  vectors  in 

a Cartesian  coordinate  system  (x,  y,  z) , in  the  directions  x,  y,  and 

z respectively.  The  corresponding  spherical  coordinates  (r, Q , 4)  with 
+ + -*■ 

unit  vectors  e , e , and  e are  related  to  x,  y,  and  z by 
r q 9 

x * r Sin  ©Cos  ♦ 
y * r Sin  9 Sin  $ 

1 • r Cos  0 

On  Bl>  the  first  term  of  Eq.  (IV-19)  gives  a velocity  vector 


8nL 


-j  (2  Cos  0*er  SinQ  e*) 


Also,  on  Bt . one  has  r * Le  and  n * e . 

la  r r 


One  thus  has 

r xln  x v) 


a Sin  0 . 

— 2 e 

2 9 

8n  l 


8irL4 


( Sin  0 Cos  0 Cos<p  i ♦ Sin  0 Cos  0 SinQ  j -Sin  0 k J 


Placing  the  above  expression  into  the  surface  integral  of  Eq.  (IV-26), 
. 1 * 

one  obtains  -j  o as  the  value  of  that  integral.  Noting  that  the  first 
integral  on  the  ri?ht-hand  sid*  of  Eq.  ( IV“26)  is  it  . *»»<» 
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/ vdR 

*L 


(IV-27) 


It  should  be  pointed  out  that  Eq.  (IV-27)  is  valid  for  any  sufficiently 

large  but  finite  value  of  L.  Although  Eq.  (IV-26)  is  independent  of 

the  sise  of  R^,  one  may  not  consider  L to  be  infinitely  large.  The 

integral  J vdR  is  in  fact  indeterminate  if  an  0.  For  example,  the 
R oo 

r-coraponent  of  this  integral  over  the  region  r>L  is,  according  to  the 

first  term  of  Eq.  (IV-19), 
a 2rt  v 

to  J J 

$ =0  0 =0 


J -p  (2Cos^0  -Sin^  0 )Sin0drd0  d 0 
r*L 


This  integral  gives 


i 

0 =o 

which  is  indeterminate. 

I • 

It  is  simple  to  show  that  the  integrals  of  the  x-  and  y-components 

• • • • 

of  v over  the  infinite  region  r >L  are  also  indeterminate.  Also, 
the  integrals  of  the  velocity  moment  r x v over  the  infinite  region 
r>  L are  indeterminate.  To  show  that  these  latter  integrals  do  not  diverge, 
terms  of  order  r ^ in  Eq.  (IV-19)  muBt  be  considered. 

Eq.  ( IV-22 ) is  valid  only  for  three-dimensional  problems.  For  two- 
dimensional  problems,  one  has,  instead  of  Eq.  (IV-22), 


J (CO80  -Cos^Q  ) 


oo 

r*L 


1 


I 


3 5 


The  two-dimensional  version  of  Eq.  (IV-26)  is  therefore 


J vdR  • 
\ 


r -*■  + r ■¥  -*• 

j(r  x<*>  )dR  - I r x(n  x v)dB 

\ ®L 


(IV-29) 


where  R^  is  a circle  of  radius  L bounded  by  B^. 

Let  the  flow  be  in  the  x~y  plane,  with  the  x-axis  selected  to  be 
in  the  °o  direction.  On  B^,  the  first  term  of  Eq.  (IV-21)  gives  a 
velocity  vector 


2UL 


■^2  (Cos©  e^.  ♦ Sin  0eg) 


where  e^  and  eg  are  unit  vectors  in  the  cylindrical  coordinate  system 
(r,  0)  given  by 


x * rCo80 
y - rSin0 


On  B. , one  has  r • Le  and  n « e . Thus  one  has 

Lt  T T 


r x(n  x v)  - - TTsr 


2HL  Sin0  e0 


- - ♦ Cos0j) 


Placing  this  expression  into  the  surface  integral  of  Eq.  (IV-29), 
one  finds  this  integral  to  be  -jci  . Equation  (IV-29)  therefore  gives, 
for  two-dimensional  problems, 


/ *vdR  - i K 

*L 


(IV-30) 


Equation  (IV-30)  is  independent  of  the  size  of  R^,  as  long  as  L 

....  + ^ 

is  finite.  The  integral  J"  vdR  is  indeterminate  if  a Q *0  because 

IW, 

_2 

the  first  term  on  the  right-hand  side  of  Eq.  (IV-20)  is  of  order  r 

Ref.  11  gives  a formula  similar  to  Eq.  (IV-30)  but  with  the  coefficient 
for  aQ,  in  equation  (IV-30),  replaced  by  1.  That  formula  is  derived 
by  considering  a finite  cylindrical  volume  centered  on  the  z-axis,  bounded 
by  two  planes  z “ 0 and  z * d,  and  using  Eq.  (IV-23).  The  derivation 
of  the  formula  neglected  the  contribution  of  the  last  integral  in  (IV-29) 

(See  the  equation  between  Eqs.  A-3  and  A-4  of  Ref.  11).  The  formula  is 
in  error  for  any  finite  cylinders.  As  it  turns  out,  however,  if  the  cylinder 
is  infinite  in  length,  then  the  formula  in  Ref.  11  is  correct, 
as  shown  in  Chapter  V. 

5.  Velocity-Moment  Integrals  in  Large  Bounded  Regions. 

In  this  Section,  a formula  for  the  velocity-moment  integral  is  derived. 
Using  the  equation 


r x v = - 


1 

•j  r U) 


■jVx(r2v) 


(IV-31) 


one  obtains,  with  the  help  of  Stokes'  Theorem 
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where  B is  Che  boundary  of  R. 

Let  R be  R^',  Che  portion  of  the  fluid  region  bounded  internally 
by  the  solid  boundaries  Bj  and  externally  by  B^.  One  has  f row  Eq.  (IV-32) 


i /*  2* 

r x vdR  - - | J r id  dR 


R ' 
Kf 


V 


^ r2(v  x n)dB  - -j  I ^r2(v  x n)dB 

\ j *3 


(IV- 


Let  R be  Rj  in  Eq.  (IV-32),  one  obtains 


- (^r2(v  x n)dB  * j (r  x v)dR  + J*r2«dR 
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Placing  Eq.  (IV-34)  into  Eq.  (IV-33)  and  noting  that  the  normal  vector 
in  these  two  equations  are  directed  opposite  to  each  other,  one  obtains 
upon  simplification 
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where  R^  include  the  solid  regions  Rj  and  the  portion  of  the  fluid 
region  R^' . 

Let  be  a spherical  region,  or  a circular  region  in  two  dimensions, 

centered  on  the  origin  with  the  radius  L.  Let  L be  sufficiently  large 
so  that  the  total  vorticity  outside  R^  is  negligible.  One  has  then 
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The  last  integral  in  Eq.  (IV-35)  is  therefore  zero  because  of  the 
principle  of  total  vorticity  conservation  discussed  in  Section  III-2. 
One  obtains  therefore,  upon  noting  r w is  negligible  outside 
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where  d is  the  integral  of  the  second  moment  of  vorticity  defined  by 


(IV-38) 


Equation  (IV-37)  relates  velocity-moment  integral  to  an  integral 
of  the  second  moment  of  vorticity.  This  equation  is  valid  for  both  two- 
and  three-dimensional  problems. 

The  fact  that  the  last  integral  in  Eq.  (IV-35)  vanishes  can  also 
be  shown  by  noting  that  if  L is  sufficiently  large,  then  a scalar  velocity 
potential  ♦ exists  on  and  is  single-valued.  The  boundary  integral 

in  Eq.  (IV-35)  can  therefore  be  expressed  as 
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Using  Stokes'  theorem,  this  integral  is  expressible  as 
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and  is  zero  because  the  integrand,  being  the  curl  of  the  gradient  of  a 
function,  is  zero. 

It  is  not  difficult  to  show  that  the  terms  of  the  asymptotic  expression 
for  v,  Eqs.  (IV-19)  and  (IV-21)  are  all  expressible  as  gradients  of  a scalar 
function.  For  example,  for  two-dimensional  flows,  the  first  term  on  the 
right-hand  side  of  Eq.  (IV-21)  and  the  next  term,  being  respectively  of 
orders  r”2  and  r are  important.  The  first  term  is  already  written 
in  the  form  of  a gradient  of  a scalar  function,  namely  -gp  $(  V lnp  'a0^* 

The  next  term  can  be  shown  to  be  expressible  as 
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V.  AERODYNAMIC  FORCE  AND  MOMENT 

In  this  Chapter,  the  theorems  and  formulas  presented  in  Chapters  III 
and  IV  are  utilized  to  derive  general  formulas  relating  the  aerodynamic 
forces  and  moments  to  the  vorticity-moment  integrals. 

1.  Aerodynamic  Force 

Consider  the  control  volume  Rj^  bounded  externally  by  ^ and  containing 
the  fluid  occupying  the  region  R^'  and  the  solid  bodies  occupying  the 
regions  Rj.  The  momentum  theorem  gives 


^ J"p  vdR  + ^*pv(v.n] 


n)dB 


(V-l) 


«L 

where  Ffc  is  the  total  force  acting  on  the  system  within  R^,  and  p 
is  either  the  fluid  density  p or  the  solid  density  P^,  depending  on 
the  particular  region  of  interest. 

On  account  of  the  asymptotic  behavior  of  v^  the  last  integral  in 
Eq.  (V-l)  is  negligible  for  sufficiently  large  values  of  L.  The  integral 
over  can  be  written  as  sums  of  integrals  as  follows: 
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The  total  force  Ffc  acting  on  the  system  in  R^  consists  of  the 
-*• 

force  F.  acting  on  the  boundary  B_  and  the  external  force  F : 

l e 
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(V-3) 


41 


The  external  force  F is  the  force  exerted  froih  outside  the  system 

e 

and  acts  on  the  solid  bodies.  The  total  force  acting  on  the  solid  bodies 
is  F + F,  where  F is  the  aerodynamic  force  exerted  by  the  fluid  on 
the  solid  bodies.  Newton's  second  law  of  motion  gives 


N d r 
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Placing  Eqs.  (V-3)  and  (V-4)  into  Eq.  (V-2)  one  obtains 
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For  sufficiently  large  L,  the  shear  stress  is  negligible  on 
and  F is  expressible  in  terms  of  p as  given  by  Eq.  (III-21).  That 

Is 

is 
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The  asymptotic  behavior  of  v shows  that  the  term  does  not  contribute 

to  the  last  integral  in  Eq.  (V-6).  The  function  f(t),  being  independent 
of  position,  also  does  not  contribute  to  Eq.  vv-6).  One  thus  has 
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Using  Eqs.  (III-20)  and  (IV-19),  one  obtains  for  three-dimensional 


flows 
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Expressing  the  lest  integral  in  a spherical  coordinate  system  with  ci£ 

•*  > 

directed  in  the  k direction,  one  obtains 
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Using  Eqs.  (IV-27)  and  the  above  equation,  one  obtains  from  Eq.  (V-7) 
the  following  expression  for  aerodynamic  force  in  three-dimensional  flows. 


t 1 dao  Ed 

' 2*  dt  + P j-1  dt  J Vj 
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For  two-dimensional  flows,  using  Eqs.  (III-20)  and  (IV-21),  one  obtains 


N 


.dR 


(V-8) 


I 


Moment  of  Aerodynamic  Force 


The  theorems  of  moment  of  momentum  gives 
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where  lt  is  a position  vector  describing  the  line  of  action  of  the  total 


force  Ft< 


The  last  integral  in  Eq.  (V-ll)  is  negligible  for  sufficiently  large 

L.  The  total  moment  consists  of  the  moment  1L  x FL  acting  on  the  boundary 

+ + ... 

B and  the  externally  applied  moment  1 x F . The  shear  stress  is  negligible 

- 

on  B^.  Consequently  the  moment  1^  x F^  is  negligible.  One  thus  has 
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The  total  moment  acting  on  the  solid  bodies  consists  of  the  externally 
applied  moment  and  the  moment  exerted  by  the  fluid  on  the  solid  bodies. 

This  total  momentum  is  equal  to  the  time  rate  of  change  of  angular  momentum 
of  the  solid  bodies.  This  rate  of  change  of  angular  momentum  is  given 
by  the  last  term  of  Eq.  (V-12).  One  therefore  obtains  the  following  expressions 
for  the  first  moment  of  aerodynamic  force  in  both  two-and  three-dimensional 


ri  j 


problems,  by  using  Eq.  (IV-37) 
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where  1 is  a position  vector  describing  the  line  of  action  of  the  aerodynamic 
force  F . 

Equations  (V-5)  and  (V-10) express  the  aerodynamic  force  exerted  by 
a fluid  on  solid  bodies  immersed  in  the  fluid  as  integrals  of  the  first 
moment  of  the  vorticity  vector  plus  an  inertia  term.  Equation  (V-13) 
expresses  the  first  moment  of  the  aerodynamic  force  as  an  integral  of 
a second  moment  of  the  vorticity  vector  plus  a moment  of  inertia  term. 

Higher  moments  of  the  aerodynamic  force  can  be  related  to  integrals  of 
higher  moments  of  the  vorticity  vector. 


3.  Summary . 

The  preceeding  results  form  a theory  for  aerodynamic  forces  and  moments 
r which  encompasses  much  of  the  previous  aerodynamic  theories.  As  discussed 

in  Chapter  II  of  this  report,  this  theory  deals  with  unsteady  motions 
of  a fluid  which  is  at  rest  at  certain  initial  instant  of  time.  A steady 
flow,  if  it  exists,  is  considered  to  be  approached  asymptotically  at  large 

> 

time  levels.  It  has  been  shown  that  under  this  circumstance  the  effective 
vortical  regions  surrounding  and  trailing  the  solid  bodies  are  of  finite 
extent  at  any  finite  time  level.  The  theory  comprises  of  the  following 
| three  statements : 

a.  The  combined  total  vorticity  of  the  fluid  and  the  solid  bodies 
is  zero. 

b.  The  aerodynamic  force  acting  on  the  solid  bodies  is  expressible 
as  the  time  derivative  of  an  integral  of  the  first  moment  of  the  vorticity 
plus  a inertia  term. 

i 

♦» 
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c.  The  moment  of  aerodynamic  force  acting  on  the  solid  bodies  is 
expressible  in  terms  of  the  time  derivative  of  an  integral  of  a second 
moment  of  the  vorticity  plus  an  moment  of  inertia  term. 

These  three  statements  are  expressible  mathematically  as  follows: 

a 

JiUdR-O  (111-5),  (III- 

where  R*  is  the  infinite  unlimited  region  jointly  occupied  by  the  fluid  and 
all  the  solid  bodies. 
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where  d is  the  dimensionality  of  the  problem. 
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